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Abstract 

Motivated by the possibility of combining spintronics with molecular structures, we investigate 
the conditions for the appearance of spin-polarization in low-dimensional tubular systems by con- 
tacting them to a magnetic substrate. We derive a set of general expressions describing the charge 
transfer between the tube and the substrate and the relative energy costs. The mean-field solution 
of the general expressions provides an insightful formula for the induced spin-polarization. Using 
a tight-binding model for the electronic structure we are able to estimate the magnitude and the 
stability of the induced moment. This indicates that a significant magnetic moment in carbon 
nanotubes can be observed. 
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I. INTRODUCTION 



Over the last decade there has been an explosive increase of activity in two key areas of 
material science: spintronics and molecular electronics. Spintronics is based on the use of 
the spin degrees of freedom as well as the electronic charge for a number of applications Q| . 
The field has expanded significantly since the discovery of the giant magnetoresistance effect 
in magnetic multilayers 2] and has been the main driving force leading to the development 
of the present generation of magnetic storage devices. 

Although the field has already demonstrated part of its potential, it is worth noting that 
most of the proposed applications simply translate well known concepts of conventional elec- 
tronics into spin systems. The typical devices are made with molecular beam epitaxy growth, 
and lithographic techniques; a bottom-up approach to spintronics devices has scarcely been 
explored. In this respect molecular electronics provides the opposite approach P]. Here the 
basic idea is to use molecular systems for electronic applications and conventional electronic 
devices such as transistors [J], negative differential resistors and rectifiers have already 
been produced at the molecular level. 

A few experiments have attempted to combine spintronics with molecular devices. In their 
pioneering experiment 0] Tsukagoshi and coworkers demonstrated that the I-V curve of a 
carbon nanotube sandwiched between two Co contacts presents hysteresis when a magnetic 
field is applied. Such spin-valve behavior indicates spin injection into the nanotube with 
a spin diffusion length (the average distance that an electron travels before flipping its 
spin direction) of the order of 100 nm. This makes carbon nanotubes very attractive for 
spintronics applications. Other carbon structures are capable of accommodating net spin 
polarization and Coey et al have shown evidence for a strong induced magnetic polarization 
at room temperature in a graphite system with embedded ferromagnetic nanoclusters ^. 

From these experiments it emerges that, on the one hand spins can propagate in car- 
bon materials almost without flipping their direction, and on the other that the proximity 
with magnetic materials can induce spin polarization in graphite-based systems. Although 
more controlled experiments on synthetic nanostructures are needed, we believe that the 
implementation of spin physics in carbon systems is possible and it will be crucial for the 
development of smaller and more sophisticated magneto-transport devices. 

Motivated by the idea of combining spintronics with molecular structures, we investigate 
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the conditions for which an induced spin polarization appears in a low- dimensional tubular 
molecule contacted to a magnetic material. Although carbon nanotubes are the immedi- 
ate motivation for this work, our formalism is rather general and can be applied to any 
cylindrical-like structures. 

This paper is organized as follows. In the next section, we derive a general expression 
for the charge transfer that occurs when a tubular molecule is side-contacted to a metallic 
substrate. A complementary expression for the contact-induced total energy change is also 
presented. This set of expressions determines not only how the charge is redistributed when 
the tubular molecule is contacted to the substrate but also provides information about the 
stability of the transfer process. The expressions derived in section IH] do not present any 
explicit spin dependence. In section IIIII we generalize them to include the spin asymmetry 
of magnetic substrates and we demonstrate that a net magnetic moment can be induced in 
the tubular molecule. An estimation of the magnitude of the induced magnetic moment and 
of its stability is then given. 

II. CONTACT-INDUCED CHARGE TRANSFER AND ENERGY GAIN 

In order to investigate how a magnetic contact affects the spin polarization of a nanotube, 
we start by calculating the change in the electronic structure of a tubular molecule side- 
contacted to a substrate. We model the contact by introducing an electronic coupling 
between the tube and the substrate that accounts for the possible charge transfer between 
the two materials. The interaction is assumed to be only between the two lines of atoms, 
one on the tube and one on the substrate, that are in closest proximity. Spin-dependent 
charge transfer is likely to arise due to the spin asymmetry of the magnetic substrate, thus 
leading to a net induced spin polarization. 

The spin-dependent density of states is the relevant quantity to calculate and provides 
the necessary information about whether or not the contact leads to a net polarization of the 
nanostructure. We start by making no assumptions regarding the specific models describing 
the electronic structure of the system. In this way, we express the spin-polarized density 
of states in terms of single-particle Green functions matrix elements that can be calculated 
by different techniques based on model Hamiltonians. Such a model-independent treatment 
emphasizes the generality of our results and leads to a set of closed-form expressions that 
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provide a general method to investigate charge transfer between contacting materials. 

We consider infinitely long tubes of diameter d, which can be thought of as 2-dimensional 
finite-width stripes wrapped around in cylindrical shape. We also assume that N atoms are 
placed along the circumference. The tubular system shows translational symmetry along the 
longitudinal direction. Therefore, the electronic states along this direction are well described 
by a reciprocal-space wave vector that runs within the 1-dimensional Brillouin zone. Since 
the translational symmetry is broken by the line of atoms contacting the substrate, it is 
convenient to use real space coordinates along the circumferential direction. In this way, 
electronic states are identified by the pair of indexes (j, k), where k corresponds to the wave 
vector along the longitudinal direction and j labels the lines of atoms on the tube surface. 
Since the substrate also has translational symmetry along the axial direction of the tube, 
the electronic states of the substrate can be labeled by the same pair of indexes. 

In terms of the single-particle Green function, the total density of states p{E) is written 

as 



where Gjj{E, k) is the Green function of an electron with energy E moving on line j with 
wave vector k. The sum over j accounts for all the atomic lines of the tube and the substrate. 

It is convenient to define another Green function associated with the disconnected system, 
that is, the isolated substrate and tube. In this case, the translational symmetry along the 
circumferential direction is reestablished and the electronic states on the tube are usually 
described by a set of two wave vectors; the longitudinal component running continuously 
over the 1-dimensional Brillouin zone and a finite set of quantized wave vectors induced 
by the cylindrical boundary conditions. Likewise, the in-plane translational symmetry is 
also restored for the substrate. However, to be consistent with the notation in terms of 
the indexes {j,k), we label the Green function of the disconnected system as Qjji{E,k) 
describing electronic propagation between lines j and j'. To distinguish between atomic 
lines on the tube and on the substrate, we label the former by integers j < and the latter 
by J > 1. It is clear that Qjji vanishes if j and / refer to lines on different subsystems. 

The effect of the contact on the total density of states can be calculated by summing up 
the corresponding change in the Green function over all possible lines, i.e., J2j ^Gjj, where 
^^3,3 ~ ^ 3,3 ~ ^3,3 gives the variation of the Green function evaluated at line j. Consistently 
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with our notation, we define lines j = and j = 1 as those atomic hnes connecting the tube 
with the substrate, respectively. Therefore, the matrix elements of AGjj are 

^G,,j = gj,o t (1 - 6^1,1 go,o t)-^ Gi,! Goj (2) 

and 

= Gj^i (1 - 6^0,0 1 Gi,! Gofi t Gij , (3) 

for lines on the tube (j < 0) and on the substrate (j > 0), respectively. In the equations 
above, we have introduced the parameter t describing the coupling between the tube and the 
substrate. This quantity plays the role of a tight-binding-like energy-independent electronic 
hopping between the relevant overlapping orbitals on either side. The Green functions above 
are energy- and fc-dependent matrices, whose indexes may refer to orbital and spin degrees 
of freedom. 

We define Ap{E) = p{E) — Pq{E) as the density of states change between the disconnected 
(po) and the side-contacted system (p). According to Eq.(^ the variation of the total density 
of states is given by 

Ap{E) = {--)ImTr Y.T.^G,,{E,k) , (4) 

^ k j 

where the trace accounts for possible internal degrees of freedom such as spin and orbital 
indexes. By combining the cyclic property of the trace with the definition of Green functions, 
we can write Ap{E) as 

Ap{E) = -{-)imTr Y,{i - g,,,tg,,rt^Y^[^tg,,,t^ + g.^.t'^t^] , (5) 

TT ^ anj an/ 

or in a further simplified form 

Ap{E) = -(-) ImJ2^ lndet(l - 6^0,0 1 t^) • (6) 

Equation (jHI writes Ap{E) in terms of the Green function matrix elements for the discon- 
nected system {g) plus the coupling parameters t and between the tube and substrate. 
More specifically, it only depends on the diagonal elements ^0,0 and ^11, namely those where 
the connection takes place. Eq.® is therefore a convenient expression to calculate the effect 
of the coupling since it provides the required change in the density of states without the need 
of evaluating the electronic structure for the connected and disconnected system separately. 
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This is equivalent to the Lloyd's formula j9[ describing the variation of density of states due 
to a diagonal perturbation, a very useful method to treat substitutional impurities. This 
method has also been used in the study of magnetic coupling between impurities in metallic 
systems |10j. The fundamental difference in the case presented here is that the perturbation 
has an off-diagonal form representing the contact between the two structures. It is worth 
highlighting that the expression above involves no approximation and is exact for arbitrary 
values of the parameters t and t^. 

From the variation in the total density of states, we can derive the change in the number 
of electrons AA^ (at zero temperature). This is obtained by integrating Ap up to the Fermi 
level Ep, i.e., 

AN{Ef) = -(-) ^/m In det(l - ^^o.oI^f) t t^) • (7) 

k 

In the equation above the relevant matrix elements are evaluated at the Fermi level. Since 
Ap{E) is the variation of the total density of states, its integral gives the change of number 
of particles in the closed system. This is of course a conserved quantity and the equation 
AN{Ef) = allows us to calculate the value of the Fermi level. 

Another quantity derivable from the change in the density of states is the effect of the 
contact on the total electronic energy. This is a fundamental quantity whose value determines 
whether or not the perturbation in the electronic structure is energetically favorable. It is 
defined as 

AE = [^^ dEEAp{E) . (8) 
From the expression for Ap{E) in Eq.® we have that 

AE = {-) V ['^^ dEIm In det{l- Go otGiit^) . (9) 

TT ^ J-oo 

If AE < 0, the changes in the electronic structure predicted by the equations (0), © and 
are possible when the energy gain is sufficient to overcome the energy costs involved in 
the transition. 

The changes Ap, AA^ and AE describe the effect of the contact on the density of states, 
number of particles, and total energy of the entire system, i.e., the tube and the substrate. 
In order to investigate possible contact-induced spin-polarizations, one must look at similar 
changes on the separate parts. In other words, instead of evaluating the total density of 
states summed over all possible sites in the structure, we must distinguish between the 
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changes in the tube and in the substrate. Bearing in mind that global charge neutrality 
is imposed by Eq.©, any modification in the total number of particles on the tube must 
be compensated by the corresponding change on the substrate. Therefore, to calculate the 
charge transfer between the tube and substrate it is sufficient to evaluate the variation of 
number of particles on either part. We choose to focus on the tube and calculate the change 
of density of states summed over all atomic lines of the tube. Analogously to the derivation 
presented above, the change in the density of states (pt) on the tube is given by 

Ap,{E) = (1) Jm Tr - g^^^ ^0,0 t)-'^t ^1,1 . (10) 

The charge transfer to the tube (AA^f) is the integral of the above expression and is written 

as 

AN,{E) = {-)ImTr dE {I - g,^,t'^ g.^otV '^tg.^it^ . (11) 

TT ^ J~oo ah, 

Equations Q and (fTT|) are the fundamental results of this work. The first tells us whether 
or not the charge transfer is energetically favorable, and the second the amount of charge 
exchanged between the tube and the substrate. These form a closed system of equations 
written in terms of the Green functions of the disconnected system and the coupling pa- 
rameters. It is important to stress that although the complexity involved in evaluating the 
expressions for ANt and AE depends on the choice of the Hamiltonian used to describe 
the electronic structure of the system, the validity of equations and (fTTj) does not. This 
means that our expressions can be equally used with simple model Hamiltonians or with a 
full realistic description of the electronic structure. 

Equation (jlOj) for Apt is not as concise as its counterpart Eq.ijH} but it can be further 
simplified by expanding it to second order in t. This approximation is valid in the limit of 
weak coupling, which is satisfied in the case of carbon nanotubes sitting on top of transition 
metals. In fact, recent density functional theory calculations of graphite on a [001] cobalt 
surface suggests a value for the coupling parameter of t = VT/SO, where W is the width of 
the graphite 7r-band Furthermore, when the sum over k is eliminated and Apt{E) is 
integrated up to the Fermi level, ANtlEp) becomes 

ANtiEp) =Tr dE 

where po{E) and pi{E) are the density of states on the tube and on the substrate, respec- 
tively. Vm{E) = tRe[gm,m]t^ (^ = 0, 1) plays the role of an energy-dependent electronic 
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PoiE) 



dVi{E) 
dE 



Pi{E) 



dVojE) 
dE 



(12) 



potential. Vi(-E') is the potential felt by the tube due to the substrate and Vo(-E') is the 
analogous potential felt by the substrate and produced by the tube. It is worth noting 
that the potentials Vo(-E') and Vi{E) depend on the real part of the Green functions ^0,0 
and Qi^i, respectively. These two quantities are directly obtainable from electronic structure 
calculations for the disconnected system. A similar second-order expansion of Eq.Q also 
simplifies the formula for the energy change AE, which now reads as 

AE = Tr dE [po{E)V,{E) + p^{E)Vo{E)] . (13) 

J —00 

Although Eqs. ()12j) and ()13|) represent a more concise version of their respective counter- 
parts, Eqs. (fTT|) and Q, they are still in integral forms. These can be further simplified by 
replacing Vm{E) with its mean value {Vm{E)). This approximation gives rise to the following 
two expressions: 

ANtiEp) = TripoiEp) [Vi{Ef) - {Vi)] + pi{Ef) [Vo{Ef) - {Vo)] } (14) 

and 

AE = Tr{No{EF){V,) + N,{EF){Vo)} . (15) 

In this form, the expressions for the charge transfer and the respective energy gain are 
written in terms of the density of states for both the tube and the substrate, and the 
potentials Vo{Ep) and Vi{Ep), all evaluated at the Fermi level. They also depend on the 
total number of electrons Nq and A^i. 

In calculating the average (Vq) ((Vi)), the integration limits are not in the range [ — 00, Ep^ 
as in Eqs. (jl2j) and (jl3|) but start from the bottom of the electronic band pi (po)- The upper 
integration limit is common to both cases and is given by the Fermi level Ep- It is clear 
from Eq. p4|) that the sign of ANt{Ep) is fully determined by the potentials Vq and Vi. The 
side-contacted nanotube will then be electron- (hole-) doped for positive (negative) values 
of ANt{Ep). This means that the balance between (Vm) and Vm{Ep) determines the type of 
charge transfer between the structures. We have checked the results obtained by the mean- 
field equations against those predicted by Eqs. (|TT|l and and we find both a qualitative 
agreement and values of the same order of magnitude. 
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III. MAGNETIC SUBSTRATES 



The expressions presented in the previous section display no exphcit dependence of the 
electronic structure on the spin degree of freedom. However, when a magnetic substrate is 
considered the spin symmetry is broken. In this case the expressions derived in the previous 
section are still valid since an explicit spin-dependence can be added to both the Green 
functions and the coupling parameter, without loss of generality. Although general non- 
collinear spin Hamiltonians can be considered, we restrict our analysis to collinear spin in 
the two-spin fluid model. Within this model all the quantities are diagonal in the spin 
subspace and the only variation with the spin-degenerate case is that the expressions for 
the charge transfer and the energy gain are different for the two spin sub-bands. Therefore 
any induced magnetization M on the tube results from the spin-imbalance of the charge 
transfer, M = {ANJ ~ ANl)fiB, where AA^^"^ is the charge transfer for a spin a and /i^ is the 
Bohr magneton. When the substrate is magnetic the charge transfer for the majority-spin 
sub-band is different from that of the minority, leading to a net induced magnetic moment 
on the tube. 

It is worth recalling that equations (|12|) and (jl3|) are complementary and that charge will 
be transferred only if the corresponding energy gain is sufficient to outweigh the energy costs. 
In the two-spin fluid model we have to calculate the energetics of the charge transfer process 
for each spin direction. Only when the energy gain is favorable for both spins does the 

1" I 

quantity (AA^ — AA/)/iB describe the induced moment. In other words, if the energy gain 
for one spin-direction is not sufficient to surpass the energy costs, the corresponding charge 
transfer will not take place and the difference M = (AA/ — AA/)/iB will be meaningless. 
Assuming that both transitions are energetically favorable we can make use of equation (fT^ 
to write the induced moment M as 

M = /.BTr J^dE \^Po{E)^[V;\E)-V,'{E)] + [pl{E)-pi{E)]^^^ , (16) 

where the spin-polarization a is now explicitly included in the quantities describing the 
substrate. Bearing in mind that the spin-bands are split by the exchange integral A and 
neglecting possible hybridization effects, we can approximately correlate the majority- and 
minority-spin bands by Vi{E) = Vi (E — A) and pi{E) = p\{E — A). A further simplification 
can be made by expanding the substrate quantities in powers of A. In this case, the induced 
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magnetization becomes 



M = -iiB^ Tr dE \p^{E) 



d'V^iE) ^dp[{E) dVo{E) 



(17) 



dE"^ ' dE dE ^ 

Both equations ()16p and (fT7|) give the induced magnetic moment in terms of quantities that 
are directly obtainable from electronic structure calculations and provide valuable expres- 
sions to determine the contact- induced spin polarizations. Whereas the latter is valid for 
magnetic substrates whose spin bands are not significantly split, the former gives a gen- 
eral expression for the induced moment with no limitations about the electronic structure 
parameters. 




FIG. 1: Density of states po and pi and potentials Vq and Vi for a typical case: a (3,3) carbon 
nanotube (solid line) attached to a semi-infinite cubic slab (dashed line). The electronic structure 
parameters for the slab leads to a band that is 5 eV wide and centered at ei = O.757. The Fermi 
level is fixed at Ep = (vertical line). 



In order to test whether or not the magnetic contact can induce spin-polarization on the 
tube we must determine the order of magnitude of the charge transfer for a given spin band. 
This can be done within a simplified model that contains the fundamental features of the 
electronic structure of both the nanotube and the substrate. The single-band tight-binding 
model is known to reproduce well the band structure of both graphite and nanotubes of 
somewhat large diameters. The electronic hopping within the tube is described by the 
parameter 7 = 2.5 eV and is hereafter used as our energy unit. Likewise, the d-band 
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of magnetic transition metals can be described within the same model by an appropriate 
choice of band width and total number of electrons. Figure 1 shows the density of states po 
and pi as well as the potentials Vq and Vi for a typical case. We have chosen an armchair 
nanotube with = 12 atoms per ring. The substrate is modeled by a semi-infinite cubic slab 
whose electronic structure parameters lead to a band that is 5 eV wide (typical of transition 
metals) and that is centered at an arbitrary position ei = O.757. The Fermi level is fixed at 
Ep — and is represented in the figure by a vertical hue. For this choice of parameters, the 
calculated results are a charge transfer of AA^^ = — 2 x 10^ electrons//im with the respective 
energy gain of AE — —6.5 x 10~^ eV per unit cell. It is worth recaUing that the present 
calculations are for zero temperature, although the finite- temperature regime can in principle 
be accounted for by including the Fermi functions in the integrands of the expressions above. 
The negative sign of ANt indicates that for this particular band alignment, electrons flow 
from the nanotube onto the substrate. 

The specific values of ANt and AE depend on the particular alignment of the electronic 
bands of the tube and the substrate. We investigate different possibilities by changing the 
on-site potential ei of the substrate atoms. This corresponds to shifting the band center 
along the nanotube energy spectrum. In figure 2 we show the changes ANt and AE as a 
function of ei for two different tube diameters. Since the Fermi energy is kept at Ep = 0, a 
shift in the electronic band of the substrate also affects the total number of electrons A^i on 
the substrate. The figure shows that the charge transfer ANt can change sign, indicating 
that the tube can be doped either with electrons or with holes. However, despite differences 
in sign, the magnitude of the charge transfer | AA^t | does not change substantially and reaches 
values up to 8 X 10^ electrons//im. The fact that the charge transfer between the tube and the 
substrate can change sign depending on the band alignment, has important consequences 
on the induced magnetic moment. In fact, if the band splitting of the substrate is such 
that the charge transfer for the majority-spin sub-band has opposite sign to that of the 
minority-spin band, the spin balance on the nanotube is not only broken but maximized. In 
this case electrons of opposite spins fiow in opposite directions (for instance, majority spins 
will fiow from the substrate into the tube, and minority from the tube into the substrate). 
The picture shows that this is the case when the Fermi level lies close to opposite edges of 
the ferromagnetic spin bands. Half metals seem to satisfy this requirement and therefore 
are potential candidates for inducing large spin imbalance in nanotubes. In other words. 
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FIG. 2: Charge transfer (AA^^) and the respective energy gain per unit ceh {AE) for different band 
aUgnments. The parameter ei corresponds to the center of the substrate band. Sohd and dashed 
hnes refer to (3,3) and (8,8) armchair nanotubes, respectively. 

magnetic substrates made of half metals are predicted to be the best materials to induce 
a magnetic moment in a nanotube. Considering the results of figure 2 as a reference, we 
estimate that the maximum value of induced magnetization is M = 10^^ /i^ per unit cell, a 
magnitude that is experimentally detectable. The stability of this induced moment can also 
be estimated by the energy gain AE, which is in the order of lO^^eV^, as shown in figure 
2. Induced moments in tubes of increasingly large diameters are less stable and AE must 
saturate toward the value associated with a nanotube deposited on a graphite substrate. 



Furthermore, we note that the lattice relaxation [1^ induced by the charge transfer as well as 
the charging energy due to the low capacitance of carbon nanotubes must be considered 
when calculating the total energy costs of the electronic transition. Although these two 
factors are influential on whether the transfer process becomes energetically favorable, they 
are at least one order of magnitude below the energy gain obtained by the change in the 
electronic structure. In fact, the typical value for the quantum capacitance (per unit 
cell) of a nanotube is C = 0.35 e/V. For the parameters used in Figure [H this leads to a 
charging energy Ec = 0.0025 e\^, which is still substantially smaller that the energy gain of 
AE = 0.065 calculated here. 

Finally, we wish to briefly comment on the similarities between our contact-induced 
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spin-polarization effect and tfie problem of spin-injection from metallic systems. Whereas 
the latter is a non-equilibrium transport effect, the former is the result of charge and spin 
rearrangement toward the equilibrium configuration between the magnetic and non-magnetic 
materials. In other words, while in the spin-injection problem the electrons must travel long 
distances to be probed by a detector, this is not the case for the contact-induced spin- 
polarization. The two phenomena can indeed be addressed by a common formalism but 
a complete comparison would require a reformulation of our method in terms of transport 
quantities, which is beyond the scope of the present paper. 

IV. CONCLUSIONS 

In summary, we have presented a model that describes the charge transfer of a carbon 
nanotube in contact to a substrate. Closed-form expressions in terms of Green functions give 
the charge transfer and the respective energy gain associated with the transition. We have 
subsequently shown that when the substrate is magnetic the spin imbalance of the surface 
may lead to an induced spin polarization in the nanotube. Within a simple model that 
reproduces the basic features of the electronic structures of both nanotubes and transition 
metal surfaces, we were able to estimate the magnitude and the stability of the induced 
moment. This indicates that a measurable magnetic moment can be induced in carbon 
nanotubes when contacted to magnetic substrates. Finally, we have demonstrated that half 
metals are the best candidates for inducing a sizable magnetic moment in carbon nanotubes. 
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